A Wavelet Based Suboptimal Kalman Filter for Assimilation 
of Stratospheric Chemcial Tracer Observations 
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A suboptimal Kalman filter system which evolves error covariances in terms of a trun- 
cated set of wavelet coefficients has been developed for the assimilation of chemical tracer 
observations of CH A . The truncation is carried out in such a way that the resolution of 
the error covariance is reduced only in the zonal direction, where gradients are smaller. 
Assimilation experiments which last 24 days, and used different degrees of truncation were 
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but it gives explicit information on the location and scale of the most important covari- 
ance structures. This allows for adaptive schemes to truncate the wavelet expansion to a 
small number of significant terms, and potentially could provide insight into the relationship 
between important features in the error covariance and background fields. 

Recent work on Kalman filtering of chemical tracer observations (Menard et al . , 2000 and 
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The details of the Kalman filter Assimilation scheme are given in Menard et al. (2000) , and 
we give a summary here. The evolution of a chemical constituent field along an isentropic 
surface is governed by 



THE TRANSPORT ADVECTION SCHEME AND CONSTITUENT ASSIMILATION SYSTEMS 



and is solved by an operator splitting approach. The system is solved off-line using 
wind analyses V produced by the Goddard Earth Observing Data Assimilation System 
(GEOS-DAS). The discretized equations have the form 


where ^ is the mixing ratio at time £& and M& the discretized transport model from £*, to 
^k+ 1- Model error in (3) come from discretization and the neglected right hand side of 1. 
e k assumed to have zero bias and normally distributed, with associated covariance Q ? . 

Menard et al. (2000) showed that the model error is state dependent, using the form 



where C (i,j) is the same as in (5). 






where the Kalman gain matrix K*, is given by 


Kfc + i (Hfc + iP{ +1 ) T (Hfc +1 P^ +1 Hf +1 + R£ +1 ) 1 (12) 

The error covariance dynamics, (9), are computed by Menard et al. using a two-step 
scheme that corrects for variance loss. The forecast error covariance is initially computed 



0 (*) = £<**( 2 * ~ k ) 



WAVELET REPRESENTATION OF ERROR COVARIANCES 


B ^ 

<0 X! 

f— i I -i 


o 


point makes up a single column in the covariance, as shown in Fig. 2. The 
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OO 



Description of the Assimilation Scheme in Wavelet Space 


We consider two schemes for carrying out the propagation of error covariances in wavelet 
space. In the first, the analysis covariance and the transport model are projected onto a 
wavelet wavelet basis and the covariance is propagated in wavelet space. The covariance 
is then transformed back to physical space and the analysis covariance and Kalman gain 
are calculated. The second approach leaves the analysis covariance in wavelet space and 
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Project the model dynamics onto wavelet basis 



ASSIMILATION RESULTS 



are the Kalman gain and observation operator mixed physical/wavelet space. Notice that 
R £+i is in physical space, because K and K T are only right and left multiplied by the 
wavelet transform matrix respectively. 



4112 x 4112. Truncation is carried out on the finest scales first. The finest scale is always 
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Constituent assimilation was carried out for 24 days (the same initial state as in the 
previous section) using the approximation scheme in (32-38), with L=1024, 512 and 256. 
These are again compared with the benchmark assimilation which evolves the full error 



ASSIMILATION RESULTS 



wavelet coefficients respectively. Both the qualitative shape and range of values of changed 
significantly here. These more severe truncations remove the (j=9,10) scales from the error 
covariance, which correspond to resolutions of 20-40° in the zonal direction. It is likely that 
this level of truncation removes too much of both zonal direction structure and propagator 
dynamics. 


resolution in the zonal direction begins to affect the assimilated field. 
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Figure 2. Schematic representation of how the error covariance for a single point (below) 
becomes a single column of the whole error covariance matrix (above). 
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